Semiclassical ordering in the large- N pyrochlore antiferromagnet 
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We study the semiclassical limit of the Sp(N) generalization of the pyrochlore lattice Heisenberg antiferro- 
magnet by expanding about the TV — > oo saddlepoint in powers of a generalized inverse spin. To leading order, 
we write down an effective Hamiltonian as a series in loops on the lattice. Using this as a formula for calculating 
the energy of any classical ground state, we perform Monte Carlo simulations and find a unique collinear ground 
state. This state is not a ground state of linear spin-wave theory, and can therefore not be a physical (TV = 1) 
semiclassical ground state. 
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Geometrically frustrated antiferromagnets have at- 
tracted interest because their large classical ground state de- 
generacy can allow a rich variety of correlated states, in- 
cluding (at T = 0) quantum spin liquids or complex or- 
dered states. The simplest examples are nearest-neighbor, 
exchange-coupled antiferromagnets in which spins from tri- 
angles or tetrahedra that share corners 0: the kagome, 
"checkerboard", and SCGO (SrCrgpGa^-gpOig) lattices 
in two dimensions, plus the garnet and pyrochlore lattices in 
three dimensions: the pyrochlore, in particular, consists of 
tetrahedra whose centers form a diamond lattice. The Hamil- 
tonian is H = Jij^i ■ Sj, where Jy = 1 for nearest neigh- 
bors (ij) . In fact, additional terms - dipole interactions and 
anisotropies (as in Gd2Ti2 0>), magnetoelastic couplings (as 
in ZnVgOu and Z11CT2O4)- decide the order in most real ma- 
terials lilyj]. Still, the case with pure Heisenberg exchange 
is worth understanding since (i) most simulations are done for 
this case; (ii) the more realistic systems emerge from it by 
the addition of perturbations; (iii) this has motivated exper- 
imentalists to search for model systems in which the afore- 
mentioned perturbations are small; (iv) quantum effects can 
be studied without being overshadowed by classical effects. 

What is the ground state for large spin length 5? In un- 
frustrated antiferromagnets, it is just the classical ground state 
dressed with zero-point fluctuations of harmonic spin waves, 
and in frustrated cases the spin-wave zero-point energy may 
lift the degeneracy of classical ground states 0. In the py- 
rochlore case, though, a large degeneracy remains ; its res- 
olution by higher-order (anharmonic) terms in the semiclassi- 
cal (1/S) expansion requires arduous approximations 0. 

An established alternative to the spin-wave approach is to 
generalize the Heisenbergspins [with SU(2) = Sp(l) symme- 
try] to Sp(TV) symmetry |7|]: here TV is the number of flavors 
of Schwinger bosons whose bilinear form represents a gener- 
alized spin 0, with length k = 2S. The resulting mean-field 
theory (valid in the /V — > oo limit) is popular as an analytic 
approach to the 5 = 1/2 limit, since the small- k limit cap- 
tures various disordered and exotic ground states . The 
large- TV mean-field theory is also useful at large-n for this pa- 
per's problem, since it gives a simple analytical prescription 
for ground state selection: unlike the spin-wave expansion, 
here all degeneracies are (typically) broken at the lowest or- 



der [0(1/ k)] quantum correction 

However, on highly frustrated lattices this approach has the 
complication of a macroscopic (exponential) number of de- 
generate saddle-points not related by symmetry, so it is un- 
known a priori which of these should be expanded around; 
this was handled till now by limiting the investigation to order- 
ing patterns of high symmetry and small magnetic cells, or by 
enumerating all saddle -points in a small finite system I9i ll(ll . 

In this letter, we develop an effective Hamiltonian ap- 
proach to this question. The pertinent saddle-points are la- 
beled by arrangements of valence bond variables, and we ob- 
tain a simple formula for the large-TV mean-field energy of any 
classical ground state, as a function of these variables. The ef- 
fective Hamiltonian is constructed as an analytical real-space 
expansion of loops made of valence bonds. This allows us to 
systematically search for a collinear pyrochlore ground state, 
using Monte Carlo annealing, on quite large system sizes. 
However, we also find that the pyrochlore ground state does 
not agree with even the lowest-order term in the spin-wave 
expansion, and therefore cannot give the right answer for the 
physical (TV = 1) ground state, in the large-S" limit, demon- 
strating a limitation of the large-TV approach for this case. 

Large TV mean field theory. — We begin by discussing the 
mean-field Hamiltonian derived from the Sp(TV) generaliza- 
tion of H. For the TV = 1 case we can write the spin inter- 
action in terms of Schwinger boson operators as 5, • Sj = 



bl a .bi cr 'bj cr ,bj <T , where a sum over repeated indices a and a' 
(that take values | , J.) is implied. The Hilbert space of the 
spin model is obtained by constraining the number of bosons 
on each site b\ a bi a = 2S. We can rewrite the interaction in 

t i,t 



terms of valence bonds created by the operator e aa ib\ a b 



where e 



jcr'' 



-eii = 1. An arbitrary singlet state can be writ- 
ten in terms of some arrangement of these bonds with at most 
2S bonds emanating from any lattice site. Generalizing these 
bond operators to TV-flavors allows us to put a large number 
of bonds on a link. Since the Hamiltonian acting on a state 
changes at most two bonds per link, the relative change in the 
number of bonds goes like 1 /TVS'. In the large-TV limit, their 
fluctuations are quenched. Therefore, we factorize the inter- 

'bi<j^mbjcr',rri)IN-, 

as classi- 



action in terms of valence bonds Q y - = (e a 
where the flavor index m = 1,2. 



., N. We treat Q y - 



cal quantities, to obtain the mean-field Hamiltonian 

Hmf = \J2 H^l 2 + (^' b l, m bl', m Qij + H, 



(1) 



Here Lagrange multipliers Ai have been introduced to enforce 
the constraint on boson number Nk at every site i, defining 
the generalized spin length k = 2S. In what follows, we shall 
take Xj to be spatially uniform A; = A. We are interested 
in large enough values of k to condense a flavor mode of the 
itinerant bosons, (bi a _ m ) = vNSi^mXi^, for long-range order 
to develop. The mean-field ground state energy (per flavor) 
is obtained by diagonalizing Q by a canonical Bogoliubov 
transformation: 



N 



i 

I [Tr \J\ 2 l - QtQ - N S X] 



(2a) 
(2b) 



Here N s is the number of lattice sites, and ( I2b> is the zero- 
point energy contribution of the bosons. The exact mean-field 
ground state is obtained by a constrained minimization of the 
above expression. It can be systematically approached as an 
expansion in powers of 1/k. The leading contribution to the 
energy (of order k 2 ) comes from terms in d2al , whose min- 
imization simply relates the valence bonds to the condensate 
configuration in the classical ground state(s) of the Heisenberg 
Hamiltonian H with spin size k/2. We will denote this con- 
figuration of bond variables with a superscript c: {Q^ }. The 
quantum correction (of order k) is provided by terms in ( I2bt 
for these bond configurations. 

The ground states of the classical Hamiltonian d2al consist 
of all spin configurations in which the spin vectors sum to zero 
in every tetrahedron. On general grounds we expect quantum 
corrections to select collinear ground states from the classical 
manifold 151 El [Till . We therefore restrict our attention to such 
states, in which each spin can be denoted by an Ising vari- 
able rji <G {±1}- Collinearity implies that, up to an arbitrary 
gauge transformation, Q?. = k(i}i — r)j)/2 and thus the bond 
variables are ±« for every satisfied, antiferromagnetic (AFM) 
bond, and zero otherwise. Also, A c = An for all pyrochlore 
lattice classical ground states. 

Loop expansion and effective Hamiltonian. — Next,we re- 
cast the first quantum correction to the mean field energy, 
Eq. J2bl >. for a given collinear classical ground state, into an 
effective Hamiltonian form where only some of the degrees 
of freedom remain |j. Eq. OB 

can formally be Taylor- 
expanded 
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FIG. 1 : (a)-(d) Schematic diagram of terms contributing to the con- 
stant term in the effective energy, due to TrQ 2 (a), and TrQ 4 (b,c,d). 
These are (001) projections, where the crossed squares are projected 
tetrahedra, and AFM bonds are shown in dark. All paths that do not 
contain loops, e.g. (a,b,c), can be viewed as paths on a coordination 
4 Bethe lattice, (e)-(f): Examples of the two types of paths that we 
need to count, in order to calculate the effective Hamiltonian coeffi- 
cients, as shown on a (001) slice of the pyrochlore lattice. The dashed 
lines represent bonds that connect to adjacent slices, (e) A decorated 
Bethe lattice path of length 14 contributing to F(14) and (f) A path 
of length 22 containing a loop of length 8, contributing to G(8, 14). 



Since \Qij/n\ = 1 for AFM bonds, and zero otherwise, 
Tr(QtQ/ft 2 ) m is equal to the number of closed paths of 
length 2m, composed of AFM bonds. All terms in Eq. Q 
depend solely on the structure of the network formed by AFM 
bonds. Note that since this network is bipartite, each nonzero 
element of QtQ is k 2 . 

In any collinear classical ground state, each tetrahedron has 
two up spins and two down spins, and four AFM bonds form- 
ing a closed loop (see Figs.^-d). This means that the local 
connectivity of the AFM network is identical for all states, and 
many closed paths only contribute state-independent terms to 
Eq. 0. For example, TrQ^Q = AN s k 2 , for any classical 
ground state since the only paths of length 2 involve going 
to and fro on the same bond, and each site has four neigh- 
bors which have the opposite spin (see Fig. ^). Similarly 
Tr(QtQ) 2 = (16 + 12+A)N s k 4 , where the three terms corre- 
spond to the paths shown in Figs.^,^,^1, respectively. All 
paths that do not involve loops, (e.g. those in Figs.^,^),^) 
can be viewed as paths on a Bethe lattice of coordination 4, 
and would contribute a constant term to the energy for all 
collinear classical ground states. The same is true for paths 
involving only trivial loops, in addition to the Bethe lattice 
path, as in Fig.^1. Here, a "trivial" loop is the loop of length 
4 that exists within any tetrahedron. The lowest order terms 
in expansion that contribute a state-dependant term in the 
effective Hamiltonian are for 2m = 6, since the shortest non- 
trivial loops are hexagons. 

This leads us to parameterize the effective Hamiltonian in 
terms of the various non-trivial AFM loops. 
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where {K21} are numerical coefficients, and V21 is the number 
of non-trivial AFM loops of length 21, per site. 
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coefficient analytical 



numerically fitted 



K 


-0.59684 


-0.59687 


K 6 


-3.482 xlO" 3 


-3.522x10 


K s 


-3.44xl(T 4 


-3.76x10" 


Ki 


-3.59 xlO" 5 


-4.5 xlO" 5 


K 12 


-3.8xl0 -6 


-5.5xl0" 6 



TABLE I: Coefficient values for Eq. 0, obtained analytically, and 
by an independent numerical fit to the energies in Fig. [2] 



To evaluate the coefficients {K 21}, we need to calculate two 
types of terms: (i) The number F{2m) of closed paths of total 
length 2m on a decorated (with trivial 4-loops) coordination-4 
Bethe lattice (Fig.QJ). (ii) The number G(2l, 2m) of closed 
paths of length 2(m + 1), involving a particular loop of length 
21 with decorated Bethe lattice paths emanating from each site 
along the loop (Fig.Q)- Calculating these terms is a matter of 
tedious but tractable combinatorics. We find that the functions 
F(2m), G(2l, 2m) decay rapidly with m, allowing us to sum 
them in order to evaluate the coefficients to any accuracy in 
Eq. (0}, using 



K = £ F ( 2m ) 

rn—0 



K 2 i = J2 G ( 2/ > 2m ) • ( 5 ) 



771=0 



We show the first five coefficients in Tab. H] Thus we have 
obtained an effective Hamiltonian that is parameterized solely 
by the number of AFM loops of various sizes. Note that the 
coefficients decay rapidly K21+2/K21 ~ 1/10, which leads us 
to expect short loops to be the dominant terms in the expan- 
sion. This allows us, in principle, to calculate the energy, to 
any accuracy, for any member of an infinite ensemble of clas- 
sical ground states. This represents a significant improvement 
over previous calculations that were always limited to small 
system sizes ll^t FTc^l . 

Although we derived the effective Hamiltonian for collinear 
states, it turns out that, in fact, the classical tetrahedron zero 
sum rule implies that Eq. (|4}, with the coefficients in Tab. H] 
is valid for any non-collinear classical ground state, as well, 
with the generalized loop variables expressed as sums over 
non-trivial loops 



v * = ^x E MQl li2 Qi 

(h—iu) 



'22!- 



^Qhih) ■ ( 6 ) 



Unlike the collinear case, where the elements of Q^Q could 
only take the values or k 2 , and thus each loop would con- 
tribute or 1 to the sum l|6}, in the general case, the matrix 
elements of Q^Q are complex. 

Numerical results. — To verify the validity of the effective 
Hamiltonian @, we calculated the energy for a large num- 
ber of collinear classical ground states, as well as linear spin- 
wave ground states, obtained by a random flipping algorithm 
described elsewhere We find that the energies are re- 

markably well described by , even when we cut the ex 
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FIG. 2: Calculated energies E q of 50 sample classical ground states 
(open symbols), 16 of which are harmonic spin-wave ground states 
(squares), along with Eq S , with 21 < 6 (dashed line) and 21 < 8 (solid 
line). The constant term Ko was numerically fitted (see main text). 
The inset shows the linear spin-wave energy for the same states. Al- 
though the spin-wave energy tends to be lower for states with lower 
E q , the large- N ground state need not be a spin- wave ground state. 



coefficient values of Tab. Q] but had to adjust the constant term 
Kq separately for each choice of cutoff, in order to get a good 
fit 11211 . In practice, this means that the effective Hamilto- 
nian @ is extremely useful for comparing energies of various 
states, even with a small cutoff, but requires many terms in 
order to accurately determine the energy. An independent 5- 
parameter numerical fit, to Eq. 10}, up to 21 = 12, gives the 
values shown in the right-hand column of Tab.Q] 

Now that we have an approximate formula for E q , for 
any collinear classical ground state, we can systematically 
search these states, with large magnetic unit cells, to find a 
ground state. We conducted Monte Carlo simulations using a 
Metropolis loop flipping algorithm and the effective energy of 
Eq. @, for various orthorhombic unit cells of sizes ranging 
from 128 to 3456 sites, with periodic boundary conditions. 
We find a minimum energy of E q /(Nn/2) = -0.60077A^ S 
for a family of nearly degenerate states. They are composed 
of layers, that can each be in one of four arrangements, re- 
sulting in ~ e cL states, where L is the system size, and c is 
a constant. Each of these states has Vq = N s /3, which is 
the maximum value that we find (but is not unique to these 
states), and "Ps = 23N s /6. Upon closer investigation, how- 
ever, we find that a unique ground state (depicted in Fig. [3J| 
is selected. The energy difference to nearby states is of order 
lO" 7 ^, corresponding to the 21 = 16 term. 

Discussion. — It was noted by one of us j5| (see also JT3I1 ) 
that, in the pyrochlore, the degeneracy of ground states of the 
spin-wave quantum Hamiltonian, at the lowest order in 1/S, is 
associated with a gauge-like symmetry. This symmetry char- 
acterizes the degenerate sub-manifold of collinear spin ground 
states by the condition 



11* 



-1 



(7) 



pansion @ off at 21 = 8, as shown in Fig. [2] We used the for all non-trivial hexagons. Since the spin-wave theory is 
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FIG. 3: The ground state of our large- AT theory, as viewed in a (001) 
projection. Here, light (dark) bonds represent AFM (ferromagnetic) 
bonds (unlike in Fig.0. The shown pattern is repeated along x and 
y directions, as well as in adjacent z slices. This state has a 48 site 
magnetic unit cell. 
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board lattice, the energy is lowest for states that have the 
most AFM (square) non-trivial plaquettes. Thus, the non- 
degenerate ground state is clearly the (7r,7r) state in which 
all plaquettes are AFM ficll . In the kagome case, all clas- 
sical ground states are non-collinear. However, if we limit 
ourselves to coplanar arrangements, we find that Qij has the 
same absolute value for all of the lattice bonds, but the signs 
differ depending on the chirality of the triangle to which the 
bond (i, j) belongs. Therefore, the effective Hamiltonian @, 
with the generalized variables l|6), prefers classical ground 
states with negative product of triangle chiralities around all 
hexagons. One can thus conclude that the ground state is the 
^/i x state, as large- N calculations have indeed found OJ. 
Let us also remark that our method can be generalized to long- 
range Heisenberg interactions which are relevant in the con- 
text of real materials like Tb 2 Ti 2 7 OH 

Finally, it has been suggested that the disordered (small-n) 
limit of the large-A^ approximation for the pyrochlore lattice 
also has a massive multiplicity of saddle-points lll4ll : an ef- 
fective Hamiltonian similar to this paper's could organize the 
handling of this family. 

UH and CLH acknowledge support from NSF grant DMR- 
0240953, and PS acknowledges support from the Packard 
foundation and the U.S. Dept. of Energy, under Contract No. 
W-31-109-ENG-38. 



FIG. 4: Per-site large- N energies E q calculated for various har- 
monic spin-wave ground states (bottom), compared to the per-site 
spin-wave energy obtained from an anharmonic calculation for S = 
1500 (top). In both cases the lowest energy is for a state that maxi- 
mizes the number of AFM hexagons. 

expected to be exact in the limit of infinite S, the physical 
ground state must satisfy Eq. Q. The state depicted in Fig. [3] 
however, does not. Looking at the inset in Fig. [2] we find 
that states with negative hexagon products tend to have lower 
large-iV energy than other states, since they tend to have more 
AFM loops, but this is not a strict rule. Thus it would seem 
that the N — > oo, large-K ground state cannot be the physical 
(N = 1) large-.? ground state. 

Nevertheless, if we restrict the large-A^ calculation to har- 
monic spin-wave ground states only, we find that the ordering 
of energies for various states is similar to preliminary anhar- 
monic spin-wave results @], and does predict the same ground 
state. As shown in Fig. |4j in both cases, the lowest energy 
among harmonic spin-wave ground state belongs to a state 
with the most AFM hexagons. 

The effective Hamiltonian approach that we have outlined 
here can easily be applied to other lattices. In the checker- 
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